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We perform a microscopic theoretical study of the generic properties of competing magnetic phases 
in iron pnictides. As a function of electron filling and temperature, the magnetic stripe (single-Q) 
order forms a dome, but competing non-collinear and non-uniform double-Q phases exist at the 
foot of the dome in agreement with recent experiments. We compute and compare the electronic 
properties of the different magnetic phases, investigate the role of competing superconductivity, 
and show how disorder may stabilize double-Q order. Superconductivity is shown to compete more 
strongly with double-Q magnetic phases, which can lead to re-entrance of the C 2 (single-Q) order 
in agreement with recent thermal expansion measurements on K-doped Ba-122 crystals. 
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In correlated materials in general, and unconventional 
superconductors in particular, a microscopic understand¬ 
ing of the magnetism is of paramount importance. Gen¬ 
erally, this is because a proper description of the relevant 
exchange mechanism in these materials is intimately tied 
to their basic electronic properties. More specifically, it is 
additionally shown within a wide class of models that the 
nature of the magnetic fluctuations may be closely linked 
to the emergence of the superconducting condensate.^ilt^ 
Focussing on the iron-based superconductors, the 
prevalent magnetic structure consists of collinear mag¬ 
netic stripe (MS) order with in-plane moments oriented 
antiferromagnetic (ferromagnetic) along the a (b) axis of 
the orthorhombic Fe lattice as shown in Fig. I^a). Thus, 
this configuration of moments singles out the Qi = (tt, 0) 
ordering vector (IQ), i.e. M(r’) = Mi exp(*Qi • r). An 
obvious question, however, is why the system does not 
take advantage of the enhanced susceptibility at both 
Qi and Q 2 = (0, tt) to form other magnetic phases, e.g. 
double-Q (2Q) phases consisting of superpositions of or¬ 
dering at Qi and Q 2 with M(r) = 2 exp(jQ;-r). 

This question has been studied theoretically mainly us¬ 
ing various effective field theories restricted to the vicin¬ 
ity of the magnetic transition temperature These 

works have identified two competing magnetic structures 
of the 2Q type: 1) an orthomagnetic (OM) non-collinear 
phase with nearest neighbor moments at right angles as 
shown in Fig. Sb). and 2 ) a collinear non-uniform spin 
and charge ordered (SGO) phase as shown in Fig. [^c). 
The favorable magnetic order depends delicately on the 
band structure, doping level, and interactions.®*^ 

Experimentally, the dominating magnetic order in the 
iron pnictides is the IQ MS state. This phase lowers 
the C 4 symmetry of the high-T tetragonal phase to or¬ 
thorhombic C 2 , and causes an associated splitting of the 
crystal Bragg peaks due to magneto-elastic coupling. Re¬ 
cently, several experiments have, however, reported the 
discovery of magnetic order without an associated struc¬ 
tural splitting, i.e. in the tetragonal phase,which has 
been taken as indirect evidence for a magnetic driven 
structural transition in the case of IQ MS order.E^lln the 


case of Ba(Fei_a;Mna;) 2 As 2 ,^ however, additional experi¬ 
ments have shown that Mn induce local regions of mag¬ 
netic (tt, tt) order and a crystal structure consistent with 
intertwined short-range clusters of both tetragonal and 
orthorhombic structure.l^lli^ The collective outcomes of 
these experiments have recently been shown to be cap¬ 
tured within a microscopic disorder scenario.^ 

The study of 2Q order in Ba(Fei_a;Mna;) 2 As 2 should be 
contrasted to other pnictides with out-of-plane dopants 
where a disorder scenario seems less relevant. This in¬ 
cludes for example Ba-122 doped with Na or K where 
experiments have found evidence for a phase transition 
into a long-range or dered m agnetic phase with tetrago¬ 
nal crystal struct ure.l^EsEMUl Por the case of hole doped 
Bai_a;Naa;Fe 2 As 2 this magnetic phase exists at the foot of 
the magnetic dome in the phase diagram.l^ More recently, 
Bdhmer et at, used thermal expansion measurements to 
find a tetragonal low-T phase consistent with a magnetic 
2Q phase, and additionally revealed a superconductivity 
induced re-entrant orthorhombic phase at even lower t\^ 
Collectively, these experimental findings define the fol¬ 
lowing main challenges for a theoretical description: 1 ) 
the existence of 2Q phases in restricted (intermediate) 
doping regimes limited by MS order and superconduc¬ 
tivity (SC), 2) The 2Q phases exist in a limited (inter¬ 
mediate) T range, 3) SC competes with magnetic order 
causing a lowering of Tc upon entering the magnetic 2Q 
phase from the paramagnet, and 4) SC competes more 
strongly with 2Q phases than the C 2 MS as seen by an 
upward T). jump when transitioning from the coexistence 
phase of SC and 2Q magnetism to a coexistence phase of 
MS order and SC.I^ 

Here we perform a study of the stability of, and phase 
transitions between, the competing magnetic phases ob¬ 
tained from a five-band Hamiltonian relevant to the iron 
pnictides. The Coulomb interaction is treated within 
unrestricted self-consistent Hartree-Fock, i.e. all charge 
and spin densities are allowed to vary at each sepa¬ 
rate site. These calculations constitute a comprehen¬ 
sive microscopic study of the 2Q magnetic phases, also 
in the presence of SC and disorder, and allows to ac- 
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cess the entir e T regime contrary to previous theoretical 
studiespHsHH xhe model provides an explanation of all 
four challenges outlined above. 

The starting Hamiltonian consists of a five-orbital 
tight-binding band relevant to the pnictide^^ and the 
usual multiband Hubbard-Hund interaction term. When 
mean-held decoupled this leads to the following model 


F 





ay 
av iiv 



where creates an electron at site i with spin a 

in orbital state p-. and are functions of 

the interaction parameters U, J (J = U/4), and 
the helds {cl^^Cj,^a') (see Supplementary Information 
(SI)). We diagonalize Eq. 0 unrestricted on N^: x 
Ny lattices, and self-consistently calculate the spin 
^'(^i) = and charge density 

nin) = > ’^liere l = x,z, and 

extract the ordering components Mi and ni of M(r) = 
X);M;exp(IQ; • r) and n(r) = n; exp(iQ; • r). We 
have compared the free energy of the magnetic states 
T = — TS to verify the stability of the results. 

Figure [^d,e) display two representative the phase di¬ 
agrams as obtained from self-consistently diagonalizing 
Eq. Q as a function of T and electron hlling for U = 0.85 
eV (d) and U = 0.95 eV (e). Uopt = 5.91 is de- 
hned as the electron hlling with the highest T)v as de¬ 
duced by the paramagnetic susceptibility, and the to¬ 
tal hlling is (n) = Uopt — x. The fact that the opti¬ 
mal doping level for the magnetic order is offset from 
n = 6 for DFT-generated bands is well known,!^ and 
not important for the conclusions of this paper. As seen 
from both cases, the 2Q phases, OM and SCO, exist 
at the foot of the MS dome on the electron and hole- 
doped side respectively, and whereas the transition be¬ 
tween the MS and OM phases is sharp, a more grad¬ 
ual transition takes place between the MS and the SCO 
phases as indicated by the green intermediate regions. 
Interestingly, both recent Mossbauer spectroscopy stud¬ 
ies of Sri_ 2 ;Naa;Fe 2 As 2 !^ and moun spin rotation mea¬ 
surements on Bai_a;K 2 ;Fe 2 AsJ^ found that indeed the 
2Q magnetic phase of the hole-doped system seems to be 
the SCO phase, in agreement with the phase diagrams in 
Fig.Sd, e). The colinear spin structure of the SCO phase 
was also recently verihed by spin polarized neutrons, ad¬ 
ditionally finding that the moments are oriented along 
the c-axis.li^The grey area denoted IC in Fig.[^e) repre¬ 
sents an incommensurate magnetic phase where the or¬ 
dering vectors Q 1 /Q 2 no longer faithfully represent the 
magnetic ground state of the system. The IC phase is 
absent in Fig. [^d) because the lower U leads to a van¬ 
ishing magnetization at significantly lower doping levels 
compared to Fig. [^e). 

In the remainder of this paper we focus on the U = 0.95 
eV case, and use the x = —0.09 and x = 0.02 electron and 
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FIG. 1. (Color online) Spin and charge order of the (a) IQ 
MS (M 2 = 0), (b) 2Q OM (Ml T M 2 ), and (c) SCO (Mi = 
M2). The black circles in (c) represent the Qi -|- Q2 = (rr, tt) 
charge order. (d,e) Magnetic phase diagrams as obtained from 
Eq. 0 as a function of T and filling Uopt — x for (d) U = 0.85 
eV and (e) U = 0.95 eV. The green (purple) area indicates 
regions of coexisting MS and SCO (OM). 


hole fillings, respectively, to discuss the transition from 
the MS state to the OM and SCO phases upon lowering 
T. Figures [^a,b) show the T dependence of the SDW 
components Mi = |Mi| and M 2 = IM 2 I for both Filings. 
At Tjv, Ml gradually increases while M 2 remains zero, 
signalling a second order transition into the MS state. 
In the X = —0.09 case, [Fig. 0a)], at Ti < the IQ- 
2Q transition takes place and the moments re-orient to 
form the OM state with Mi T M 2 and M 2 = Mi. As 
seen from Fig. 0a), the sudden jump of M 2 at Ti is 
compensated by a reduc tion in Mi , leaving the average 
magnetization Mj. = ^ nearly unchanged. In 

Fig.0 c) we display the T dependence of the entropy 
S{T). The small discontinuity in S at Ti [see inset of 
Fig.0 c)] agrees with a weak first order transition. 

The IQ-2Q transition taking place at a: = 0.02 is shown 
in Fig. 0b, d). As seen, in this case the second component 
M 2 continuously increases below Tq < Tjv in a second 
order fashion. The increase of M 2 is again compensated 
by a decrease in Mi. In this case, however, M 2 is aligned 
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FIG. 2. (Color online) (a) T evolution at a: = —0.09 
of the magnetic components Mi (t), M 2 (a) and Mr = 
l-y/Mj + M| (o), and (c) the entropy S(T). The color 
changes in Mr and 5 represent the magnetic phase transi¬ 
tions shown also in Fig.[^e). (b,d) The same as (a,c) but for 
X = 0.02 with SCO order. (e,f) summarize the T-dependence 
of the magnetic structure and the expected associated lattice 
symmetry. 


with Ml, and the spin order remains collinear. As soon 
as M 2 > 0, a small charge order is induced at Q 3 = 
Qi + Q 2 = ( 7 r, 7 r) which scales with Mi • M 2 and thus 
increases gradually (see SI for more information). In the 
range Ti < T < Tq, where M 2 is increasing, the system 
is still C 2 symmetric (0 < M 2 < Mi) but has developed 
characteristics of the SCO state, such as the (tt, tt) charge 
order [see Fig. [^c)]. Finally at T < Ti the pristine Ci 
symmetric SCO phase with M 2 = Mi is formed, but we 
find that at the lowest T < Tq a re-entrance to a weakly 
C 2 symmetric 2Q phase occurs (0 < M 2 < Mi). As 
can be seen from S{T) in Fig.j^d) weak thermodynamic 
signatures are expected throughout the T range. The T 
evolution of both transitions and their associated lattice 
symmetries are summarized in Figs. [le,f). 

Next we compare the electronic properties of the mag¬ 
netic phases, MS, OM, and SCO, by focussing on the 
three different fillings indicated by the black stars in 
Fig. e). In order to illustrate the different nesting con¬ 
ditions, we first show in Figs, [^a)-(c) the FS in the PM 
state in the folded Brillouin zone (BZ) where the X and 
Y centered elliptical electron pockets /3i and ^ 2 , and the 
M centered 7 hole pocket all fold on top of the F point 
(see SI for further details). In the ordered state, energy 
gaps can open at the crossings of the bands connected 
by the magnetic ordering vector, weighted by the ma¬ 
trix elements a”(k) = (n|/ik) which relate orbital and 
band states. This is apparent in Fig. [^a,d) where the 
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FIG. 3. Fermi surface in the folded BZ, —7r/2 < kx,ky < 
7 r/2, in the PM state for (a) x = —0.09, (b) x = —0.04 and 
(c) X = 0.02, and in the magnetic states (d) OM (Mi = 
0 . 22 ; M 2 = 0.2*), (e) MS (Ml = 0.52; M2 = 0), and (f) 
SCO (Ml = M2 = 0.43). (g-i) Band structure along the 
momentum path X — T — Y shown in (d) for the (g) OM, (h) 
MS, and (i) SCO phase. The main orbital contributions are 
shown by purple: d^z’, green: dyz\ orange: d^y', cyan: d^2. 


weakly nested /Si and the outer hole pocket a 2 are gapped 
around kx = upon the OM formation. The rest of 
the SDW gaps are opened below the Fermi energy ep, 
and the reconstructed bands are seen from Fig. [^g) to 
contain considerable orbital mixture. As the filling is 
increased in Figj^b), multiple electron-hole band cross¬ 
ings get closer to the ep, and additional nested areas 
appear at the FS, for example the ones connected by 
the 7 and Pi pockets. This enhances the spin suscepti¬ 
bility at Q1/Q2, which naturally leads to larger SDW 
order and more pronounced energy gaps and FS recon¬ 
struction in the MS state. The resulting FS in Fig|^e) 
exhibits hole-like Dirac cones along the AFM direction 
and a hole pocket at F of mainly dxz character resulting 
from the mixing of ai and 0 : 2 - Evidently, since the MS 
state singles out the Qi ordering vector, the spectrum 
becomes C 2 symmetric as shown in Fig. [^e,h). Because 
MS breaks the dxzjdyz degeneracy an associated ferro- 
orbital ordering nxz > nyz results in a splitting of the 
bands at the F point as seen from Fig. ^h). No such 
splitting takes place in the 2Q states with Mi = M 2 (at 
least in the absence of spin-orbit coupling). Finally, as 
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FIG. 4. (a) Modified phase diagram from the region indicated 
by the dashed box in Fig.[^e) in the presence of SC order, (b) 
Magnetic anisotropy 6m = (Mi —M 2 )/(Mi + M 2 ) versus T for 
four different x with (without) SC order shown by the solid 
(open) symbols. The arrows mark Tc. (c,d) Fermi surface at 
the cross in (a) for MS (c) and SCO (d) order (without SC) 
using the same orbital color code as in Fig. 


the filling is increased further, large nested areas of ai 
and 7 with the electron pockets appear at the FS as seen 
from Figj^c). For the present band, at this filling the 
FS-nesting is the strongest of the presented cases, and 
the resulting SDW order parameter and the gaps are the 
largest. As seen from Figj^f), ai and Pi become fully 
gapped, and similarly most of the 7 pocket except small 
pieces around which hybridize with the outer¬ 

most a 2 pocket. The dependence of the band reconstruc¬ 
tion on the interaction parameters is found in SI. 

Motivated by the recent experimental discovery of the 
effects of superconductivity on the magnetic states,^ 
we have included SC order to the model by the follow¬ 
ing BCS term Ubcs = where 

Aii,.(fc) = Y.k',ap^ia{k - k'){c-k'fiiCk'aP). A recent 
Landau order parameter expansion has also been used 
to study this problem close to The effective pair¬ 
ing vertices — fc') are obtained from the RPA spin 

and charge susceptibilities in the PM state with leading 
s± symmetry (see SI for details). In order to study the 
effects of SC on both the IQ and 2Q magnetic states, 
we focus on the region of the phase diagram outlined by 
the dashed box in Fig. [^e), and self-consistently solve 
the associated Bogoliubov-de Gennes equations including 
both magnetic and SC order parameters. The modified 
phase diagram shown in Fig. Qa) exhibits a noticeable 
effect of SC on the boundaries between both magnetic 
states as seen by comparison to the dashed gray lines in¬ 
dicating the transition lines from Fig. I^e) without SC 
order. (The phase diagram corresponding to Fig. m 
including superconducting order is supplied in the SI.) 


Below Tc, where both the magnetic and SC order pa¬ 
rameters are non-zero, the MS region expands at the ex¬ 
pense of the SCO phase. This effect is explicitly shown 
in Fig. I^b) by the evolution of the magnetic anisotropy 
Sm = (Ml — M 2 )/(Ml -I- M 2 ), which is a measure of the 
C 2 symmetry breaking, i.e. Sm = 0 {Sm = 1 ) for the 
SCO (MS) state, and Sm > 0 for mixed C 2 states with 
0 < M 2 < Ml. As seen from Fig. |^b), without SC Sm 
gradually evolves from Sm = ^ Sii high T to = 0 at 
low T, with a transition that sharpens with increasing 
X. In the presence of SC order, however, Sm is pushed 
up, as indicated by the shaded regions in Fig. |^b), and 
the magnetic order is driven towards the MS phase. This 
effect is particularly pronounced in the regions of large 
A^i/(fc) closer to the MS phase, and is consistent with 
recent experiments.^ 

The Tc line shown by the filled black dots in Fig. |^a) 
evidently exhibits a clear drop across the MS-SCO tran¬ 
sition. This reduction of is directly caused by the 
emergence of the SCO state as verified by the signif¬ 
icantly higher Tc (empty circles) found by a separate 
calculation with the magnetic order forced to the MS 
type. The stronger competition between SC and mag¬ 
netic order can be explained by a reduced density of 
states N{e) at ep in the SCO phase (compared to MS), 
N 2 q{^f) 0.64AfiQ(ei;’). In addition, the dominant SC 

pairing is the intra-orbital d^y element caused by strong 
7 -/ 3 i FS nesting (see SI). However, as seen explicitly 
from Figs, [^d) the FS in the SCO phase (without SC) 
contains significantly less dxy orbital character (orange 
points) compared to the corresponding MS FS shown in 
Fig. B c). In summary Tc is reduced from the PM state 
into the magnetic phase {NiQ{eF) QA2NpM{f^F) and 




FIG. 5. (a,c) Magnetization at low T and x = —0.09 in the 
absence of an impurity, and (b,d) in the presence of a single 
impurity at the center (red dot) with (b) F = 5 eV, and (d) 
F = 2 eV. The associated Fourier components Mi and M 2 
are also shown (the vectors are normalized to (a) M = 0.16, 
(b) M = 0.43, (c) M = 0.2, and (d) M = 0.68.) 
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N 2 Q{f-F) 0-‘2QNpM{iF)), but enhanced from the SCO 
phase into the MS phase at lower T in agreement with 
the experimental finding of Ref. [T51 

We end with a brief study of disorder effects on the 
three magnetic states (without SC) by including the term 
Tiimp = in Eq. Q, relevant for a non¬ 

magnetic impurity at site i*. Figure summarizes our 
main finding: impurities enhance the 2Q phase by nu¬ 
cleating local SCO islands. The left column shows the 
homogeneous magnetic phase without the disorder po¬ 
tential (MS and OM), and the right column the resulting 
SDW order in the presence of the impurity. As seen, the 
charge potential acts like a seed for the charge order asso¬ 


ciated with the SCO phase, and pushes the magnetic or¬ 
der towards the collinear magnetic SCO structure. This 
can be more clearly seen from the Fourier components 
Ml and M 2 also depicted in Fig. 

In summary, we have presented a detailed microscopic 
study of competing magnetic phases in iron pnictides. 
We have mapped out the phase diagram both in the ab¬ 
sence and presence of SC, and found qualitative agree¬ 
ment with recent experimental findings. 

We acknowledge useful discussions with A. V. 
Chubukov, I. Eremin, R. M. Fernandes, J. Lorenzana, 
and financial support from a Lundbeckfond fellowship 
(grant A9318). 
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Supplemental Materials: ’’Competing magnetic double-Q phases and 
superconductivity-induced re-entrance of C2 magnetic stripe order in iron pnictides” 

Here, we provide the details of the tight-binding band and the mean-field decoupled Hamiltonian both in real and 
momentum space. We show the orbitally resolved order parameters in the single-Q and double-Q phases. In addition, 
we show the dependence of the band and the Fermi surface reconstruction on the interacting parameters. Finally, we 
provide the details of the RPA-generated couplings used to generate the superconducting order. 


I. MODEL 

The starting Hamiltonian consists of a five-orbital tight-binding band relevant to the pnictideP^^, 


(SI) 


The operators (ci/io-) create (annihilate) an electron at site i in orbital state /r with spin cr, and ^0 is the chemical 
potential which adjusts the filling. The indices /r and v denote the five iron orbitals d^z, dyz, d^y , and dz'^- 

The five orbital tight-binding band and Fermi surface for filling x = —0. is shown in Fig. Here, x refers to 
a filling of (n) = nopt — x where Uopt = 5.91 is the optimal doping for magnetism for this band. Thus x = —0.09 
corresponds to the undoped case with (n) = 6.0. 

The interacting part of the Hamiltonian is described by the multi-orbital onsite Hubbard model 


Tdint — U ^ U ^ T (ff d) 




ifi 


+ d -b J ^ipLij^ip,a^ivaCivc! ^ 






(S2) 


with U' = U — 2J, J' = J, and J = 17/4. 

We mean-held decouple Eq. (S2) for all helds {c^ipaCjua') which leads to the following mean-held Hamiltonian 


= E (it 

ijfliy 



(S3) 


where 


+ U{n,yjz) + U'{n,,^) + (U' - J){n,,,)] - {U' - 


jm. 


t) + J + J ^ 


(54) 

(55) 




FIG. SI. (a) Band structure along the high symmetry directions (F = (0,0), X = (tt, 0), and M = (tt, tt)), and (b) Fermi 
surface with main orbital character for the undoped system. 
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Eq. (S3) is diagonalized by a unitary transformation 
unrestricted fields are obtained self-consistently 


IfJ, 


and 




and the following 


n 

(S6) 

n 

(S7) 

n 

(S8) 


(S9) 


n 


for all sites i, j and orbital combinations /r, v. Here En denote the eigenvalues, and / is the Fermi function. From 
these fields we obtain the spin and charge configurations of the final solution in real space 


= E -k (cEc^/.T)) > 

(SIO) 


(Sll) 

n{r) = E ■ 

(S12) 




In order to readily study the electronic properties we also solve the above model in momentum space with the 
mean-fields + qi)), where qi — {0,Qi,Q2,Qi + Q 2 } = {^Oj dij 92 , 93 }- The mean-field Hamiltonian in 

momentum space takes the following form 


+ivr 


E 


V 


IF^ 




WY 



e'^"(fc + 9i) 



IFf^ 




+N^'' 





^’'{k + qY 

irr 











e"(fc + 93) 
+K'' 



wY 





-irr 

-W^'' 



+N^'' 


i^Yk + qi) 


-VF^ 

+iVo^" 



H.c. 

^^Yk + q2) 

-wY 


^^''{k + qY 


+<" / 


where 

c];^(fc + 9i) cl^{k + q 2 ) c\^{k + q3) c\^{k) c^^ik + qi) c\^{k + q 2 ) + 93 )) , (S14) 


and the summation is done in the reduced Brillouin zone —7r/2 < ky < tt 12. 

The entries below the main diagonal were not included for clarity but are obtained by the transpose conjugate of 
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FIG. S2. Orbitally resolved (a) magnetic OPs for d^z, dyz and d^y, and (b) charge OPs for d^z and dyz as a function of 
temperatnre for x = —0.09. 



kT[eV] 




FIG. S3. Orbitally resolved (a) magnetic OPs for d^z, dyz and d^y, (b) charge OPs for d^z and dyz, and (c) (rr, tt) charge OPs 
for all orbitals as a fnnction of temperatnre for x = 0.02. 


the upper triangular matrix. All entries of the mean-field Hamiltonian are defined below 


e''{k + qi)=e^’'{k + qi)- (S15) 

Nr = [Un^ + (2t/' - Jho] + i-U' + 2J)nr + J'n^, (S16) 

wr=- jMj') - u'ur - j'Mr, (si?) 

Nr = [C/n^ -b (217' - J)r] + {-U' + 2J)nr + (S18) 

Nr=Sy.. (-[/< - JO - u'fir - j'nr, (si9) 

wr = - u'Mr - j'Mr, (S20) 

Nr = i-Un^ - Jns) “ U'nr - J'O- (S21) 


Using the unitary transformation -I- qi) = + 9*) = En (where Z = 0,1) 2,3) the 











4 


mean-fields are then self-consistently obtained from the relations 


! 3 


fc(7 kn l—O 

^r = E("M<.(fc)c-(fc+93)) = E| E KikMAk)+u7;{k)ui{k)] + T. [<;(fcX(fc)+<;(fc)C(fc)] [/- 

k(T kn —1,4 1—2,3 J 

Mr = Y.a{clAk)£.Ak + qi)) = i2\ E Kik)uUk)-ur{k)ufM + T. K*ik>Uk)-urik)ulik)] 


k(T 


kn [ 1 — 1,2 
I 


Z=3,4 


Mr = Y.ar^rkrc.a{k-rqr) = Y.\ E + E K;(fcX(^)-<;(fc)«L(fc)] E- 

ka kn —1,3 Z—2,4 

/ 3 

< = E^(4(fc)c-(fc)) = EE - <;(fcx(fc)] /- 


k(T 


kn l—O 


^r = E^(4(fc)c-(fc+93)) = E| E [<(^x(fc)-<(fc)<w] + E K;(fc)<(^)+<;(fc)C(fc)] 

ka kn [^Z —1,4 Z—2,3 J 

^r = E(4(fc)c-(fc+9i)) = E| E Kik)uUk)+ur{k)ui{k)] + e [<(fc)<(fc)+<(fc)C(fc)] 

ka kn [^Z —1,2 Z—3,4 J 

^2^" = E(4(fc)c-(fc+92)) = E1 E Kikhiik)+ur{k)ufM + E K*ik>Uk)+urwuUk)] \ /„, 


ka 


kn I Z=l,3 


Z=2,4 


where the abbreviation /„ = f{En{k)) has been used. 


II. ORBITALLY RESOLVED ORDER PARAMETERS: Mj^ AND nf 


The Figs. [S^ and |S3| show the orbital content of the magnetic and charge order parameters versus T corresponding 
to the results presented in Fig. 2 of the main text. The largest magnetic order parameter (OP) components Mj^ 
correspond to the best nested orbitals, dxz/dyz for x = —0.09 in Fig. [S^a), and dxy for the higher Filing of x = 0.02 
in Fig. S3 ’a) where the 7 —A nesting has improved. The remaining orbital components are not shown for presentational 


simplicity. The ferro-orbital order (ng^ > nT) can be seen across the PM-MS transition in Figs. [S^b) and [S^b), 
and it collapses upon formation of both OM and SCO order. In addition, a small = Qi + Q 2 = (tt, tt) charge order 
develops as soon as M 2 > 0 in the SCO phase. Figure [S^c) shows the orbital character of this charge order, and the 
dominant contributions arise from the dxz 


and dyz orbitals. 


III. BAND RECONSTRUCTION AND INTERACTION PARAMETERS 


The band reconstruction taking place in the ordered state obviously depends on the amplitude of the magnetic OPs 
Ml and M 2 , which in turn depend on the interacting strength U. This is illustrated in Fig. |S4| where the reconstructed 
FSs in the (a) OM, (b) MS and (c) SCO states are plotted as a function of U. For comparison the paramagnetic FS 
is shown in the first panel of each case. As seen, the resulting reconstructed Fermi surface depends significantly on U. 
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(a) 



FIG. S4. Fermi surfaces for different values of the interacting parameter U in the (a) x = —0.09 and OM, (b) x = —0.04 and 
MS, and (c) x = 0.02 and SCO. 
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FIG. S5. Spatial dependence of the intra-orbital effective pairing constants r{J'^(rij) in eV from the central site, (a) (b) 

dxz, (c) dxy, (d) d^ 2 _y 2 and (e) d^ 2 . 


IV. SUPERCONDUCTING PAIRING VERTEX 


The multi-orbital pairing vertex in the singlet channeP^l is calculated from the RPA spin- ^ind charge- 

susceptibilities 


r;;(fc-fc',o) = 


U^Xs (k - k', 0)U® + -U^ - ^U’^Xc ik - k', 0)U‘= -b -U' 


(S22) 


J pq 


where and are 5x5 matrices identical to those of Ref. 


The real-space pairings are then obtained by 
Fourier transforming equation (S22|, r^J;^(rij) = X]qexp(jq • (rj — rj)). We retain all orbital combinations up 
to next-nearest neighbor sites to calculate the superconducting order parameter = J2a/3 r^a(rij)(cj/ 3 ^CiQ,^).l^The 
intra-orbital effective pairings are shown in Fig. |S5[ 


V. PHASE DIAGRAM INGLUDING SUPERCONDUGTIVITY IN THE LOW-U LIMIT 

In addition to the phase diagrams shown in the main text, we have also calculated the phase diagram in the presence 
of competing superconductivity at C/ = 0.85 eV shown in Fig. |S6| The corresponding normal state phase diagram 
for this case is shown in Fig. 1(d) of the main text. The lower value of U pushes the magnetic structure to lower 
values of the doping which prevents the occurrence of the IC magnetic phase, and the paramagnetic superconducting 
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FIG. S6. Phase diagram showing the magnetic and superconducting phases as a function of T and filling Uopt — x for U = 0.85 
eV. 


phase directly merges with the SCO C' 4 -magnetic phase in this case. Note that for the particular parameters used to 
generate Fig. |S6| there is no superconductivity-induced re-entrance of the C 2 phase, which we attribute to the pairing 
being too weak to cause a switch of the preferred magnetic structure. 
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